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I. INTRODUCTION 



At sufficiently high temperatures, the quark-confinement property of QCD is expected to be lost so that hadrons 
dissolve into a plasma of quarks and gluons. This quark-gluon plasma state must have existed in the early Universe, 
and its experimental detection is being actively pursued through relativistic heavy-ion collisions. Basic information 
that is required to explore physical phenomena in the quark-gluon plasma is its equation of state, namely the energy 
density and pressure as a function of temperature. For this reason a number of lattice QCD studies of equation of 
state have been made An important progress in this effort is the recent work of the Bielefeld group Q in which 
a systematic continuum extrapolation was carried out for the equation of state for pure SU(3) gauge theory. Using 
the standard plaqucttc action, they calculated bulk thermodynamic quantities on lattices with the temporal extent 
N t — 4, 6 and 8, and extrapolated the results to the continuum limit N t — > oo, assuming that the data at N t — 6 and 
8 follow the leading extrapolation formula. 

To understand the quark gluon plasma in the real world, this work has to be extended to full QCD with dynamical 
quarks. This is a difficult task due to significant increase of the amount of computations needed for full QCD 
simulations. One approach to lessen the computational cost is to employ improved actions designed to have reduced 
lattice cut-off effects, and hence should allow reliable continuum extrapolation from coarser lattice spacings compared 
to the case for the standard unimproved actions. Thermal properties of several improved pure gauge actions have 
already been studied At present, however, no extrapolation of thermodynamic quantities to the continuum 

limit has been made with improved actions. 

In this article we report on our study of the continuum limit of equation of state for pure SU(3) gauge theory with 
an improved action determined from an approximate renormalization-group argument . Simulations are carried out 
on I6 3 x 4 and 32 3 x 8 lattices, and the energy density and pressure are calculated by the integral method |J. We 
find that the results extrapolated to the continuum limit agree well with those obtained from the standard action in 
Ref. This provides us with a cross-check of the final results in the continuum limit, and also provide support for 
the validity of assumptions behind the extrapolation procedures. 

This paper is organized as follows. In Sectio n O , we summarize the basic formulation and our notations. Some 
details of our simulations are given in Section OS We define our choice of the temperature scale in terms of the 
string tension in Section IV, and examine scaling of the critical temperature in Section |v|. In Section [Vl] we present 



our results for equation of state obtained at Nt = 4 and 8, and their continuum extrapolation. A comparison of our 



results with those obtained from the standard action is also made. In Section VII we briefly discuss results obtained 



with the operator method H . We end with a brief conclusion in Section VIIl 



II. METHOD 

The partition function of a finite temperature SU(3) lattice gauge theory is defined by 

Z=f[dU]e 0S °, (1) 
where f3 = 6/g 2 is the bare gauge coupling, and S g the lattice gauge action. Denoting the k x I Wilson loop in the 



(^,f)-plane at a site x as W** l (x) = iRcTr 



rifcxz loop ^ 



the renormalization group (RG) improved gauge action 



we use is given by |7) 

S 9 = c J2 W™(x) + Cl £ W™{x) (2) 

with Co = I — 8ci and c\ = —0.331. On a lattice with a size N 3 x N t and the lattice spacing a, the temperature T 
and physical volume V of the system are given respectively by 

T = l/(N t a), V = (N s a) 3 . (3) 

We calculate the energy density e and pressure p using the integral method || . For a large homogeneous system, 
the pressure is related to the free energy density / through 

p=-f = ^\nZ. (4) 
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Using an identity In Z = (S g ), we then have 



P_ 

T 4 



= / dp' AS , (5) 

Pa Jf3 



where 

AS^N?{(S) T -(S)o) , (6) 

with (S)t the expectation value of the action density S — Sg/N^N at temperature T. The zero-temperature 
expectation value (S)o is introduced to subtract the vacuum contribution, which is conventionally computed on a 
symmetric lattice with the same spatial volume N% x N s . Once the pressure is known, the energy density can be 
computed using 

e ~ 3p =As, (7) 



T 4 dT 

where 



dT = ~° = M (8) 



is the QCD beta-function. 



III. SIMULATIONS 

The fundamental quantity used in the integral method is the action difference AS defined by (|^). In order to 
calculate (S)t and (S)o, we need to simulate both asymmetric (iV 3 x N t ) and symmetric (iV 3 x N s ) lattices. The 
spatial lattice size N s should be sufficiently large to suppress finite size effects. Past finite-size studies [|lo]|ll]J| suggest 
that the condition N s /N t > 3 is the minimum requirement. As with the work of the Bielefeld group for the plaquette 
action 0, we choose N s /N t — 4 and perform simulations on 16 3 x 4 and 32 3 x 8 lattices as well as on 16 4 and 32 4 
lattices for a set of values of j3 around and above the critical point. 

Gauge fields are updated by the pseudo-heat-bath algorithm with five hits, followed by four over-relaxation sweeps; 
the combination of these updates is called an iteration. We always start from the completely ordered configuration, 
and perform 20 000 to 36 000 iterations after thermalization on asymmetric lattices, and about 10 000 iterations on 
symmetric lattices. We find these number of iterations to be sufficient for a statistical accuracy of 2-3% or better for 
the value of AS. The statistics of our runs are compiled in Table |. 

We measure Wilson loops and Polyakov loop at every iteration. Errors are determined by the jack-knife method. 
The typical bin size dependence of the jack-knife error for the action density is shown in Fig. [j]. The errors are almost 
constant over a wide range of bin sizes, and we adopt the bin size of 1000 iterations for asymmetric lattices and 500 
on symmetric lattices. 

In Table |H we list the expectation value of the action density (S)t calculated on asymmetric lattices of size 16 3 x 4 
and 32 3 x 8 and that for (S)o on symmetric lattices of size 16 4 and 32 4 . 



IV. TEMPERATURE SCALE 



In order to determine the temperature T — l/(N t a(j3)), we need to compute the lattice spacing a as a function of 



the gauge coupling f3. We use the string tension of static quark potential to fix the scale of this relation. In Table III 
we compile results for the dimensionless string tension ay/a obtained with the RG- improved action PJl2]|. We fit 
these results by an ansatz proposed by Allton jL3|, 

(ay/P}(/3) =/(/?) (l + c 2 a(p) 2 + c 4 a((3) 4 + ...)/c , a(fl) = j^r , (9) 

where f(/3) is the two-loop scaling function of SU(3) gauge theory, 

--^J* 1 11 10? 

f(P = 6//) = (W) ex P (-^) , bo = ^ , h = ^ , (10) 
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and c n (n — 2,4, ■ • •) parameterize deviations from the two-loop scaling. Truncating the power corrections at n = 4 
and choosing f3\ — 2.40, we obtain from this fit, 

c = 0.524(15) , c 2 = 0.274(76) , c 4 = 0.105(36) (11) 

with x 2 /dof = 0.356 for 4 degrees of freedom. As shown in Fig.|| the fit curve reproduces the data very well. 
With this parametrization, the temperature in units of the critical temperature T c is given by 

with (3 C the critical coupling. The beta-function is obtained by differentiating the left hand side of (^|) with respect 
to the lattice spacing a, keeping a constant. 



V. CRITICAL TEMPERATURE 



We determine the critical coupling (3 c (N tl N s ) for the deconfinement transition on an x N t lattice by the peak 
location of the susceptibility x of the Z(3)-rotated Polyakov line. The values of f3 c for Nt = 3, 4, 6 have been reported 
in Ref. S. In order to compute /3 C for Nt — 8, we perform an additional simulation of 24 000 iterations at /3 = 2.710 
on a 32^x 8 lattice. The p dependence of x is calculated by the spectral density method H]. We estimate the error 
by the jack-knife method with the bin-size of 2000 iterations. The values of f3 c (N t ,N s ) for finite N s are summarized 



in the second column of Table IV 



Calculating the critical temperature requires an extrapolation of (3 c (N t , N s ) toward infinite spatial size N s — > oo 
for each N t . For the first-order transition of the pure gauge system, the spatial volume dependence of f3 c (N t ,N s ) is 
expected to follow Q 

N? 

p D (N t ,N a ) = c (,Nt,oo)-c(N t )j±. (13) 



It has been reported in Ref. || that results for Nt = 3 and 4 reasonably satisfy ( |13| ) with c(N t ) — 0.122(54) for Nt = 3 
and 0.133(63) for Nt — 4. An approximate scaling of the coefficient c(N t ) motivates us to apply ( ^3|) for N t — 6 (as 
was made in Ref. M) and also for Nt — 8, adopting the value 0.133(63) for the coefficient. Substituting values of 
(3 e (Nt, oo) in the parametrization of the string tensio n (p|) , we calculate T c j ' yfo = 1/ '(Ntay/o). We tabulate results of 
this analysis in the third and fourth column of Table 



[V 



In Fig. H we plot the results for T c j ' \fa as a function of 1/Nf (filled circles). Also shown are the values previously 
reported in Ref. j| (open circles) and those for the plaquette action from Ref. (l(| (open squares). A slight difference 
between the present results and those from Ref. || for the same action stems from the fact that an exponential 
ansatz ^/aa = Aexpf— BP) with fit parameters A and B was adopted in the previous work, which deviates from 
the parametrization (0). We think that the present parametrization gives a better estimate of er, being theoretically 
consistent with the asymptotic scaling behavior for large /3. The results for the plaquette action is obtained with a 
parametrization 

{ayfij){p) = /OS) ( 1 + 0.2731 a{l3f - 0.01545 a(/3) 4 + 0.01975 a(/3) 6 )/0.01364 (14) 

with Pi = 6.0. 

We observe that the new value of T c j y/a for the RG-improved action for Nt — 8 is consistent with the previous 
results for N t = 3, 4, 6 The difference in this ratio obtained for the two actions, however, still remains. Making 
a quadratic extrapolation in 1/N t , we find T c jyfa = 0.650(5) for the RG action, which is 3% higher than the 
value 0.630(5) for the plaquette action (Tq |. A possible origin of the discrepancy is systematic uncertainties in the 
determination of the string tension for the two actions, which differ in details. We consider that checking an agreement 
beyond a few percent accuracy, as is needed here, would require the generation and analyses of potential data over 
the relevant range of lattice spacings in a completely parallel manner for the two actions, which is beyond the scope 
of the present work. 
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VI. EQUATION OF STATE 



A. Results for RG-improved action 

Our results for AS at Nt = 4 and 8 are shown in Fig. ||. In order to integrate AS* in terms of (3 to obtain the 
pressure, we have to make an interpolation of the data points. At large /? where T > 2T C is satisfied, we fit the points 
by a perturbative ansatz, 

AS = p + p + ^ + "-' (15) 

truncating the series at the order P~ 4 . The absence of the linear term in perturbation theory can be checked easily. 
For AS at lower f3- values corresponding to T < 2T C , we perform a cubic spline fit with the requirement that the curve 
smoothly joins to the large-/3 fit curve at T = 2T C . The interpolation curves are shown in Fig.|j. 

We use these curves to evaluate the integral for pressure in (^). The lower limit of integration is chosen to be 
/3q = 2.20 (Nt = 4) and 2.63 (Nt = 8). The results are shown by solid lines in Fig. ||[ Combining these results for p 
with those for e — 3p computed using (Q), we also obtain the energy density e, which we show in Fig. |^. 

The statistical error Sp(T) for the pressure, plotted at representative points in Fig. ^, is evaluated from the contribu- 
tions Sip(T) at each simulation point pi which is estimated by the jack-knife method. Since simulations at different Pi 
are statistically independent, we compute the final error by the naive error-propagation rule, Sp(T) — yj^i ^Pi(^) 2 j 
summing up all the contributions from Pi smaller than p corresponding to the temperature T. The error Se(T) of 
the energy density is calculated by quadrature from the error of 3<5p(T) and that for e(T) — 3p(T), the latter being 
proportional to the error of AS*. 

We observe in Figs. || and || that the energy density and pressure exhibit a sizable increase between N t = 4 and 8. 
This increase is opposite to the trend for the plaquette action, but it is consistent with the prediction of the leading 
order perturbative result shown by horizontal lines at the right of the figure. The values from the integral method, 
however, overshoot those from perturbation theory toward high temperatures, particularly for Nt = 4 for which the 



perturbative value is quite small. We discuss this point further in Sec. VII. 

We now extrapolate the results for energy density and pressure to the continuum limit Nt — * oo. The RG-improved 
gauge action has lattice discretization errors of 0(a 2 ). Therefore, at a fixed temperature in physical units, we expect 
deviations of thermodynamic quantities from the continuum limit to be 0(\/Nf): 

- T = p,t. (16) 



T A ) Nt Wcont N{ 



Extrapolating the results for e and 3p at Nt — 4 and 8 with this form, we obtain the continuum predictions drawn by 
solid lines in Fig. [?]. 



B. Comparison with results for the plaquette action 

We compare our results with those of the Bielefeld group obtained with the plaquette action j^] . Care is needed in 
this comparison since they used a scheme different from ours to fix the temperature scale: their scheme is based on 
the requirement that the critical temperature T c is independent of the temporal size Nt. 

Since a difference in the scale can sizably affect results for thermodynamic quantities p7| , we first examine the 
possible influence of this difference. For this purpose we reanalyze the raw data of Ref. B] for the action density 
employing the scale parametrization (^) . The results for p and e are shown by solid lines in Fig. |^ for the temporal 
sizes Nt = 6 and 8 used by the Bielefeld group for the continuum extrapolation. Compared with their original results, 
drawn by dashed lines, the influence of the scale determination is well within the statistical error of 1-3%. 

We extrapolate the pressure and energy density obtained with the scale ( [hi] ) to the continuum limit according to 
(filf). This leads to the dash-dotted curves shown in Fig. |^. Errors are evaluated in the same way as for the case 
of the RG-improved action. We observe that the curves for the RG-improved action (solid lines) and those for the 
plaquette action (dash-dotted lines) are in good agreement, within the error of 3-4%, over the entire temperature 
interval shown. This is highly non-trivial since results for the two actions differ significantly at finite lattice spacings. 

In Fig. we compare the pressure at Nt — 4 from the two actions to the result in the continuum limit. We note 
that the continuum results for the RG-improved action are approached from below, while those for the plaquette 
action from above. We also find that the magnitude of deviation from the continuum limit is comparable for both 
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actions. The form of the RG-improved action we employed is determined so as to best approximate the renormalized 
trajectory after a few RG transformations, within those limited actions with a maximum of 6-link loops. Therefore, 
low-momentum modes, with a momenta smaller than the inverse of several lattice spacings, are improved. On the 
other hand, a momentum scale which is significant at high temperatures is T = l/N t a on finite- N t lattices. From 
Fig. ^, it appears to be required to add further terms in the action in order to reduce the cut-off effects for high- 
momentum modes with momentum > l/4a, thereby improving the behavior of the pressure at high temperatures on 
N t — 4 lattices. 



VII. COMPARISON WITH RESULTS WITH OPERATOR METHOD 



In Sec. VIA we noted that the pressure and energy density calculated by the integral method exceed the values 



corresponding to the perturbative high temperature limit. This is a puzzling result, especially for pressure; while p/T 
has to decrease at high temperatures to agree with the perturbative result, we expect it to be an increasing function 
of temperature since it is given by an integral of AS which is generally positive. The discrepancy is particularly large 
for Nt = 4 for which the leading order perturbative results on the lattice are quite small compared to those in the 
continuum as first noted in Ref . . In Table [v| we list the perturbative value of pressure on a N^ x Nt lattice in 
units of the free gluon gas value in the continuum for N t = 4-12 for the RG-improved and plaquette actions. 

In order to further examine this problem, we calculate thermodynamic quantities in an alternative way using the 
formulae of the operator method Q9| given by 
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— A = -Nt[c s (g){{S s ) (S)o)-c t (g)((S t ) (S) )], 



^ = ~-N?P(g)[(S s + S t )-2(S) }. 



T 4 3T 4 



(17) 
(18) 



where S s and St are the spatial and temporal part of the action density, and the asymmetry coefficients are defined 
by 



Cs(g) = 1- 9 



2 dg s 2 



/ \ , , 2 d 9t 2 
ct{9) = 1 + .9 



(19) 



The one- loop values of the asymmetry coefficients for the plaquette action have been long known fll9| , and preliminary 
values for the RG-improved action have recently been reported [f^of . 

We compare results for the energy density and pressure from the integral method and the operator method with 
one-loop asymmetry coefficients in Fig. [Io| for Nt — 4 and 8. For both types of actions, the values for the operator 
method lie above those for the integral method, and the difference diminishes with increasing Nt. For Nt = 4 for 
the RG-improved action, in particular, we do not observe any indication of decrease toward the perturbative high 
temperature limit, both with the integral and operator methods, at least within the temperature range where we have 
results. 

A possible source of the discrepancy is breakdown of perturbation theory due to the infrared divergence |2l| of the 
theory In the continuum there are non-perturbative contributions to free energy beyond 3-loop level. This problem 
should also exist on the lattice, and the magnitude of non-perturbative contributions may vary depending on the 
choice of lattice actions. 



VIII. CONCLUSIONS 



In this article we have presented results on the equation of state for a pure SU(3) gauge theory obtained with an 
RG-improved gauge action. The continuum result for the energy density and pressure show an agreement with the 
results of the Bielefeld group for the plaquette action within the error of 3-4%. This provides a concrete support for 
the expectation that continuum results are insensitive to the choice of lattice actions. 

We also found that the energy density and pressure for finite Nt overshoot the perturbative high temperature limit. 
Understanding the origin of this behavior shall be explored in the future. 
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lattice P ^iterations 

16 3 x 4 2.20 - 2.30 26 000 

2.32 - 3.20 20 000 

16 4 2.20 - 3.20 10 000 

32 3 x 8 2.60 - 3.80 36 000 

32 4 2.60 - 3.80 12 000 



TABLE I. Statistics of our runs. 



a 
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1 fi 3 -v A 


ifi 4 

ID 


^9 3 v 8 


o 9 4 
oz 




ii.oozyoo^i4i ) 


ll.DOzDO t \ o I ) 






9 9^n 


1 1 8^/1 1 Q9^ 


1 1 S^Rfi^f 1 3^ 
ll.oOOoOO^loO J 






9 970 




1 1 QQ1 17^ 
ll.yol40O^ll / ) 






2.300 


12.053019(249) 


12.039327( 87) 






2.320 


12.123597(183) 


12.108015(141) 






2.350 


12.221928(144) 


12.206421( 93) 






2.400 


12.373140(177) 


12.360168(129) 






2.500 


12.645009(144) 


12.636627(102) 






2.600 


12.885558(111) 


12.880320( 93) 


12.8802888(225) 


12.8802498(198) 


2.650 






12.9923472(240) 


12.9923202(252) 


2.700 






13.0990662(657) 


13.0987350(177) 


2.750 






13.2009177(255) 


13.2000711(171) 


2.775 






13.2497922(213) 


13.2489255(162) 


2.800 


13.299015(138) 


13.296786( 99) 






2.850 






13.3897845(186) 


13.3890450(186) 


2.900 






13.4780103(183) 


13.4773656(201) 


3.000 


13.644306(141) 


13.643502( 63) 


13.6438083(141) 


13.6433112(162) 


3.200 


13.938774( 87) 


13.938480( 69) 


13.9385442(174) 


13.9383060(150) 


3.400 






14.1934539(120) 


14.1933099(150) 


3.600 






14.4165000(162) 


14.4164202(114) 


3.800 






14.6135793(141) 


14.6135304(117) 



TABLE II. Expectation value of action density (S) for our runs. 



f3 a^/a lattice # of conf. Ref. 

2.1508 0.5054(93) 9 3 x 18 400 

2.2827 0.3864(32) 12 3 x 24 200 

2.40 0.3096(54) 16 3 x 32 50 

2.5157 0.2559(23) 18 3 x 36 100 

2.60 0.2313(58) 16 3 x 32 50 

2.70 0.1963(34) 16 3 x 32 100 

3.20 0.1029(19) 32 4 50 @ 
TABLE III. Results for string tension obtained with the RG improved action. 
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A" 3 x N t 


Pc(N t ,N s ) 


)3 c (JVt,oo) 




12 3 x 3 


2.1528(9) 


2.1551(12) 


0.665(10) 


16 3 x 4 


2.2863(10) 


2.2879(11) 


0.654(4) 


18 3 x 6 


2.5157(7) 


2.5206(24) 


0.654(5) 


32 3 x 8 


2.7103(32) 


2.7124(34) 


0.652(6) 



TABLE IV. Critical coupling of the deconfinement transition for the RG improved action on an TV 3 x Nt lattice, infinite 
spatial volume extrapolation and the ratio T c /^a for infinite spatial volume. Allton's parametrization of string tension is 
employed to fix the temperature scale. 



p(N t ,N s )/psB 


N t 


4 


6 


8 


10 


12 


RG 


0.1971 


0.7086 


0.8213 


0.8734 


0.9024 


plaquette Q 


1.4833 


1.1697 


1.0748 


1.0398 


1.0229 



TABLE V. Perturbative high temperature limit of pressure on the lattice in units of its continuum value (psb/T 4 = 87r 2 /45) 
Results for the case N s /Nt =4 are listed. 
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FIG. 1. Bin size dependence of the jack-knife error of the action density 
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FIG. 2. String tension as a function of gauge coupling. Solid line represents a fit to Allton's parametrization. 
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FIG. 3. Tc/y/a as a function of l/iV t 2 . Open circles are results reported in Ref. ||. Open squares are values for the plaquette 
action 
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FIG. 5. Pressure for Nt = 4 and 8. The dashed horizontal line on the top-right represents the leading order perturbative 
value in the high temperature limit in the continuum, and solid and dotted lines are the corresponding lattice values for Nt — 8 
and 4. 
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FIG. 6. Energy density for Nt = 4 and 8. Meaning of horizontal lines is the same as in Fig. 
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FIG. 7. Equation of state in the continuum limit for the RG-improved action (solid lines) and for the plaquette action 
(dash-dotted lines). The latter is obtained with the Allton's parametrization of string tension using raw data in Ref. 0|. 
Dashed horizontal line on the top-right shows the free gluon gas value in the high temperature limit. 
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FIG. 8. Energy density (lower figure) and pressure (upper figure) for the plaquette action for N t = 6 and 8. Solid lines uses 
the Allton's parametrization of a for scale and dashed lines are original results of the Bielefeld group using a different scale 
fixing scheme. 
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FIG. 9. Pressure at N t = 4 from the RG and plaquette actions. For comparison, we also plot the result in the continuum 
limit obtained with the plaquette action using our choice of the scale discussed in the text. 
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FIG. 10. Results for energy density and pressure for N t = 4 and 8 obtained with the operator method using one-loop values 
for the asymmetry coefficients as compared with those from the integral method drawn by dashed (Nt = 4) and solid (Nt = 8) 
lines. Upper figure is for the RG-improved action, and lower figure for the plaquette action. 
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